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This  report  summarizes  research  performed  in  Contract 
P4i|620-7^-C-0048  during  the  period  January  15  through 
July  15,  197^.  The  research  was  Jointly  sponsored  by^the 
Space  and  Missile  Systems  Organization  and  the  Air  Force 
Office  of  Scientific  Research  (APSC),  United  States  Air 
Force.  The  Air  Force  program  monitors  were  Captain  A.R. 
Hunt  and  Second  Lieutenant  T.  Hopkins  of  SAMSO  and  Dr. 

M.  Rogers  of  AFOSR.  Mr.  W.  Portenier  maintained  technical 
llason  with  the  Aerospace  Corporation. 

Participants  in  the  study  were  D.C.  Wilcox,  principal 
investigator,  and  T.L.  Chambers.  Manuscript  preparation 
was  supported  by  B.A.  Wilcox,  W.A.  Coonfield,  and  G.J. 
McCornock.  Professor  P.G.  Saffman  of  the  California  In- 
stitute of  Technology  also  made  important  contributions 
during  the  course  of  the  study. 

A technical  note  based  on  results  obtained  in  this  study 
has  been  submitted  for  publication  in  the  AIAA  Journal. 

The  note  is  entitiled,  "Turbulence-Model  Transition  Predic- 
tions" and  is  authored  by  D.C.  Wilcox.  Material  presented 
in  the  note  includes,  with  the  exception  of  surface  rough- 
ness effects,  all  significant  results  presented  in  Sections 
2 and  3« 

A second  AIAA  Journal  publication  is  in  preparation.  The 
paper  will  be  based  on  the  surface-roughness  analyses  of 
Sections  3 and 


ABSTRACT 

An  accurate  and  efficient  method  has  been  developed  for 
predicting  transition  from  laminar  to  turbulent  flow.  The 
method  is  based  upon  a phenomenological  turbulence  model 
originally  devised  by  Saffman  for  fully  turbulent  flows. 

The  model  is  modified  to  account  for  low  Reynolds  number 
phenomena  and  is  used  to  make  a priori  transition  predictions 
Primary  emphasis  focuses  upon  blunt-body  flows;  predictions 
are  also  made  for  constant-pressure  boundary-layer  flows  to 
provide  as  broad  a data  base  as  possible  to  assess  model 
accuracy.  The  method  yields  excellent  agreement  between 
computed  and  experimentally  measured  effects  of  freestream 
turbulence,  surface  roughness  and  suction  on  incompressible 
flat-plate-boundary-layer  transition.  Model-predicted  ef- 
fects of  surface  roughness  on  transition  in  the  vicinity 
of  the  stagnation  point  of  a sphere-cone  body  virtually 
^ duplicate  measured  effects.  The  method  also  yields  close 

agreement  between  calculated  and  measured  surface-cooling 
and  nose-radius  effects  on  sphere-cone  transition.  Based 
on  results  obtained  in  the  study,  a blunt-body  transition- 
prediction  algorithm  for  use  by  weapon-system  design  en- 
gineers has  been  devised. 
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NOTATION 
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DEFINITION 


J 


k 

M 

N(k/e) 

p.pt 


^stag 


^N 

^®k*  ^®x* 


Re,] 

Re. 


S(ku.^/v) 
T,  T^ 

T« 

u 


V 

X 


U. 


y 

a,  a* 


Skin  friction 

Specific  heat  at  constant  pressure 
Volume  coefficient 
Specific  turbulent  energy 

Index  for  planar  (J«0)  or  axisymmetric  (J*l)  flow 
Roughness  height 
Mach  number 

Surface-roughness  function 

Static  pressure,  total  pressure 

Laminar  and  turbulent  Prandtl  number 

Local  heat  transfer,  stagnation-point  heat  transfer 

Radial  body  coordinate 

Nose  radius 

Preestream  unit  Reynolds  number 

Reynolds  number  based  on  roughness  height, 
plate  length,  momentum  thickness 

Turbulent  Reynolds  number 

Neutral-stability  Reynolds  number 

Empirical  parameter  in  Saffman  turbulence  model 

Arc  length 

Surface-roughness  function 

Static  temperature,  total  temperature 

Turbulence  intensity  at  boundary-layer  edge 

Velocity  component  in  x direction 

Boundary-layer-edge  velocity,  freestream  velocity 

Velocity  component  In  y direction 

Coordinate  lying  along  a solid  body 

Coordinate  normal  to  a solid  body 

Empirical  parameters  In  Saffman  turbulence  model 
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NOTATION  (Cont.) 

DEFINITION 

Values  of  o,  a*  for  fully  turbulent  flows 
Empirical  parameters  In  Saffman  turbulence  model 
Boundary-layer  thickness,  displacement  thickness 
Kinematic  eddy  viscosity 
Momentum  thickness 
Karman's  constant 

Empirical  parameter  In  Saffman  turbulence  model 
Molecular  viscosity 
Kinematic  viscosity 

Empirical  parameter  In  Saffman  turbulence  model 
Fluid  mass  density 

Empirical  parameters  In  Saffman  turbulence  model 
Angle  from  centerline  for  spherical  body 
Turbulent  dissipation  rate 
Specific  turbulent  dissipation  rate 
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Boundary-layer  edge 
Transition  point 
Body  surface 
Freestream 
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1.0  INTRODUCTION 

An  understanding  of  boundary-layer  transition  phenomena  in 
blunt-body  flows  is  essential  to  the  design  of  reentry- 
vehicle  nosetlps.  Transition  sensitivity  to  many  complicated 
effects  must  be  defined;  these  effects  Include  surface  rough- 
ness, surface  cooling,  mass  transfer,  nosetlp  geometry,  free- 
stream  Mach  number,  and  freestream  turbulence  level.  Although 
experimental  work  ' has  attempted  to  quantify  the  importance 
of  these  various  effects,  a reliable  analytical  tool  for  pre- 
dicting complex-flow  transition-point  location  remains  to  be 
developed.  Since  flight  conditions  often  cannot  be  simulated 
accurately  in  a wind  tunnel  and  since  applicable,  retrievable 
flight-test  data  are  rare,  accurate  analytical  tools  are  needed 
to  facilitate  extrapolation  from  wind-tunnel  conditions  to 
real  flight  conditions. 

Although  analytical  tools  for  predicting  transition  from  laminar 
to  turbulent  flow  have  Improved  significantly  in  recent  years, 
transition  remains  one  of  the  least-understood  phenomena  of 
fluid  mechanics.  Presently,  analytical  studies  focus  on  the 
following  four  approaches: 

1.  Semiemplrical  Formulas 

The  simplest  approach  is  to  start  with  a physical  fact 

or  argument  and  devise  a formula  containing  adjustable 

parameters;  a classic  example  of  this  method  Is  given 

0 

by  Van  Driest  and  Blumer  . Such  formulas  can  be  very 
useful  for  parametric  studies  and  for  correlating  ex- 
perimental data  for  a given  class  of  flows,  but  a lack 
of  universality  generally  limits  th^*  utility  of  this 
approach. 

2.  First  Principles 

The  most  rigorous  tact  that  can  be  taken  is  to  seek 
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time-dependent  Navler-Stokes  solutions.  While  pre- 

9 

llmlnary  steps  have  been  taken  In  this  direction  , 
required  computing  times  are  currently  too  lengthy 
to  make  such  computations  practicable  for  engineer- 
ing applications. 

3.  Linear  Stability  Analysis 

The  approach  which  has  received  the  greatest  amount 
of  attention  Is  linear  stability  analysls^®”^^. 

While  some  Insight  Into  the  transition  phenomenon 
has  attended  this  work,  predictions  often  differ  from 
experimental  observations.  Furthermore,  linear  anal- 
ysis determines  stability  of  a flow  to  Infinitesimal 
disturbances  only  and  Is  Inapplicable  when  Initial 
flow  perturbations  are  of  finite  amplitude.  Finally, 
In  a linear  stability  analysis,  complicating  effects 
such  as  wall  roughness,  wall  colling,  mass  transfer, 
and  freestream  turbulence  level  and  scale  consider- 
ably Increase  the  approach’s  mathematical  complexity. 

^ . Phenomenological  Turbulence  Models 

A relatively  new  method  offers  many  advantages  over 
the  approaches  listed  above;  this  method  uses  the 
Reynolds-averaged  equations  of  motion  subject  to  a 
set  of  closure  hypotheses  suitable  for  accurate  com- 
putation through  transition.  This  approach,  on  the 
one  hand,  is  applicable  to  arbitrary  amplitude  dis- 
turbances and,  on  the  other  hand,  in  a simple  and 
natural  way  can  account  for  the  complicating  effects 
cited  above.  Recent  progress  with  phenomenological- 
turbulence-model  equations  indicates  this  approach  Is 
sensible,  l.e.,  that  adequate  closure  approximations 
can  Indeed  be  determined.  Using  turbulence-model 

equations  in  which  the  Reynolds  stresses  depend  upon 
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flow  history,  Donaldson  , Jones  and  Launder  and 
Wllcox^^  have  shown  that  such  equations  accurately 
predict  abrupt  transition  from  laminar  to  turbulent 
flow  for  constant-pressure  boundary  layers. 

The  turbulence-model  approach  Is  taken  In  the  present 
study,  where  the  primary  objective  has  been  to  develop  an 
accurate  method  for  predicting  transition  In  blunt-body 
flows.  Including  the  effects  of  surface  roughness,  sur- 
face cooling,  mass  transfer,  blunt-body  geometry,  free- 
stream  Mach  number,  and  freestream  turbulence.  To  accom- 
plish this  objective,  a series  of  transition-prediction 

computations  has  been  performed  based  on  the  Saffman 

19 

turbulence-model  equations  . To  provide  a broad  data 
base  for  assessing  turbulence-equation  accuracy,  the  cal- 
culations Include  transitional  flat-plate-boundary-layer 
(FPBL),  pipe,  and  channel  flow  in  addition  to  the  more- 
pertinent  blunt-body-flow  computations.  Section  2 dis- 
cusses the  model  equations  and  the  method  by  which  transi- 
tion predictions  are  made;  qualitative  features  of  the 
model-predicted  transition  mechanism  are  also  described. 
Section  3 presents  transition  predictions  for  constant- 
pressure  boundary-layer  flows  and  includes  comparisons 
with  experimental  data.  Section  ^ summarizes  transition 
predictions  for  blunt-body  flows;  based  on  the  results 
presented  In  Section  a new  engineering  correlation  Is 
devised  which  can  be  used  for  predicting  nosetlp  transi- 
tion. Results  of  the  study  are  summarized  In  Section  5. 


2.0  FORMULATION 


The  basic  analytical  approach  taken  in  the  study  for  pre- 
dicting transition  Is  presented  in  this  section.  Subsection 

2.1  presents  the  turbulence  model  equations  and  the  method 
by  which  the  equations  are  used  to  predict  transition.  Model 
revisions  required  for  Improving  transition-prediction  ac- 
curacy are  described  In  Subsection  2.2.  The  section  concludes 
with  a discussion  of  the  model-predicted  transition  mechanism. 

2.1  TURBULENCE  MODEL  EQUATIONS 

The  turbulence  model  equations  which  form  the  basis  of 

17 

the  present  study  were  originally  devissd  by  Saffman  for 

Incompressible  flows.  In  subsequent  development  efforts^®*^^, 

Wilcox  has  cast  the  model  equations  In  a form  suitable  for 

compressible  flow  applications.  The  Saffman  model  Is  a two- 

equation  model  of  turbulence;  l.e.,  two  nonlinear  diffusion 

equations  are  solved  in  addition  to  the  conservation  (mass, 

momentum,  energy)  equations  In  computing  a given  flowfield. 

20 

In  terms  of  mass-averaged  mean  quantities  , the  coupled  set 
of  equations  describing  compressible  boundary-layer  flows 
over  planar  (J*0)  or  axisymmetrlc  (J=l)  surface  are  as  follows 


Conservation  of  Mass 


•^(r^pu)  + gy(r'^pv)  = 0 


(1) 


Conservation  of  Momentum 


(2) 


Conservation  of  Energy 


pu^(CpT)  + PV^(OpT)  . + („+e/n)(|H)‘ 


3y 


Specific  Turbulent  Energy 

'=“11  + '’''If  ■ ['>*  l|fl  - * ’'ly> 

* 7T  jf  ('') 

Specific  Turbulent  Dissipation  Rate 
pu|^  + pv||-  * t^“l|yl  ” Bpfi]p«^ 

+ ?r  57  (M+«/n)|fi]  (5) 

In  Equations  (1-5),  x and  y are  orthogonal  coordinates  with 
X lying  along  the  body  and  y being  normal  to  the  surface;  r 
Is  the  radial  body  coordinate  (see  Figure  1).  Velocity 
components  In  the  x and  y directions  are  denoted  by  u and  v. 
The  quantities  p,  p,  and  T denote  pressure,  density,  and 
temperature.  Molecular  viscosity  Is  the  quantity  y while 
the  eddy  viscosity  Is  given  by  e/fi,  the  ratio  of  specific 
turbulent  energy,  e,  to  specific  dissipation  rate,f2.  Lami- 
nar and  turbulent  Prandtl  numbers,  and  Prtp,  appear  In 
Equation  (3)  while  seven  empirical  parameters  (a,a* , B , 6* ,o , 
o*,C)  are  contained  In  Equations  (4  and  5).  Subscript  e de- 
notes the  value  of  a quantity  at  the  boundary-layer  edge. 

With  the  exception  of  C,  values  of  the  empirical  parameters 
and  the  turbulent  Prandtl  number  are  regarded  as  universal 
constants  for  fully  developed  turbulent  flows,  and  their 


-5- 


values  have  been  established  by  general  arguments  based  on 
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well-documented  experimental  observations  for  such  flows 
The  value  of  C has  been  selected  by  comparison  of  model  pre- 
dictions with  experimental  data  for  a variety  of  supersonic 

l8 

and  hypersonic  flows  . In  all  past  calculations,  values 
for  these  constants  have  been  given  by: 

Pr^  » 0.89 
c « 2.50 

0.  * = 0.50 

0*  = 0.09 
a*  = a*  = 0.30 
0 = 0.50 

0.15  < 0 < 0.18 

a = = a*[0/0*  - 4oK^/a*]  J 

where  k = 0.^1  is  the  Karman  constant. 

In  solving  Equations  (1-5),  values  of  u,  T,  e and  U are 
prescribed  at  the  boundary-layer  edge.  The  no-slip  velocity 
boundary  condition  is  Imposed  at  y = 0 while  either  sur- 
face temperature  or  surface  heat  flux  is  given.  Turbulent 

19 

energy  vanishes  at  y = 0 and,  for  perfectly  smooth  surfaces 
y^n  » 2Ou/0p^  (the  boundary  condition  for  0 appropriate  on 
a rough  surface  is  discussed  in  Subsection  3.2). 

The  Saffman  turbulence  model  has  been  incorporated  in  a 
boundary-layer  program  developed  at  the  NASA  Langley  Research 
Center  ; the  modified  program  is  known  as  EDDYBL  (see  Ap- 
pendix). In  using  EDDYBL  to  make  boundary-layer  transition 
predictions,  specific  turbulent  energy  and  specific  dissipa- 
tion rate  are  held  constant  at  the  boundary-layer  edge. 
Turbulent  energy  is  set  to  zero  throughout  the  boundary  layer 
at  a point  near  the  plate  leading  edge.  Solution  of  the 
parabolic  system  is  accomplished  by  marching  in  the  stream- 
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wise  direction.  Some  entrainment  and  molecular  diffusion 
of  turbulent  energy,  e,  into  the  boundary  layer  initially 
occurs;  however  little  or  no  turbulent-energy  amplification 
occurs  for  a plate-length  Reynolds  number  below  a critical 
value  Rejj  , signifying  existence  of  laminar  flow.  Then, 
when  Re^  ^ abrupt  Increase  in  e is  observed,  fol- 

lowed by  an  asymptote  to  a value  characteristic  of  fully 
developed  turbulent  flow.  The  transitional  regime  is 
readily  identified  as  the  range  over  which  e Increases  from 
its  initially  low  level  to  its  much  higher  value  in  the 
turbulent  regime.  The  transitional  regime  can  also  be  iden- 
tified from  the  numerical  data  by  locating  abrupt  changes 
in  quantities  such  as  momentum  thickness,  shape  factor,  and 
skin  friction. 

2 . 2 LOW-REYNOLDS-NUMBER-EFFECTS 

In  an  earlier  study the  Saffman  model  showed  promise  of 
being  an  accurate  tool  for  transition  prediction.  Using 
the  values  of  the  empirical  parameters  given  in  Equation  <6), 
transition  was  predicted  for  Rayleigh  shear  flow.  Although 
experimental  data  are  not  available  for  transitional  Ray- 
leigh flow,  the  computed  transition  appeared  realistic  as 
computed  Raylelgh-flow  properties  and  measured  FFEL  proper- 
ties were  in  close  quantitative  agreement.  The  first  calcu- 
■■atlon  in  the  present  study  was  for  FF3L  flow  so  that  more- 
direct  comparison  could  be  made  between  model  predictions 
and  experimental  data. 

Figure  2 shows  computed  skin  friction,  c^.,  as  a function  of 

Re  for  an  Incompressible  FPBL.  The  freestream  value  of 
* -9  2 

e is  10  Ug,  where  is  freestream  velocity.  As  shown, 
the  predicted  transition  begins  at  Re^^  = and  ends  at 

Re  = 10®.  The  predicted  value  of  Rex^  of  ^'id*  is  much 
Icv^rer  than  the  measured  for  very  low  freestream  tur- 
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Figure  2.  Skin  friction  as  a function  of  plate— length  Reynolds  number 
for  an  incompressible  FPBL. 


bulence  levels;  Schubauer  and  Skramstad-' , for  example, 
indicate  that  should  be  nearly  3 *10®  for  e©  '''  10“'^Ug. 

Assuming  a,o*,6,8*,o,a*,C  are  independent  of  Reynolds  number 
leads  to  the  inaccuracy.  As  noted  above,  the  values  for 
the  empirical  parameters  in  Equation  (6)  are  presumed  valid 
for  fully  developed  turbulent  flows;  however,  there  is  no 

a priori  reason  why  these  parameters  should  be  independent 

” 22 

of  Reynolds  number.  In  fact,  Rotta  argues  that  such 

parameters  depend  directly  upon  the  turbulence  spectrum; 

since  the  turbulence  spectrum  will  be  quite  different  for 

low-Reynolds-number  turbulence,  Rotta 's  argument  suggests 

that  a, a*,  etc,  should  be  Reynolds-number-dependent.  Other 
15 

Investigators  have,  in  fact,  introduced  a functional  de- 
pendence of  similar  parameters  upon  turbulent  Reynolds  num.- 
ber  Re^  = e/fly.  In  the  present  study,  transition  is  found 
most  sensitive  to  the  values  of  a and  a*;  decreasing  o* 
tends  to  delay  transition  while  the  ratio  of  a to  a*  fixes 
the  width  of  the  transition  zone;  therefore  to  improve  tran- 
sition-prediction accuracy,  an  assumption  is  made  that 

a«  « a«  [1  - (1-A)(l-Re^/R^)K(1-Re^/R^)]  (7) 

a/a*  = a^/a*  (8) 

where  H(x)  is  the  Heaviside  steofunction  and  and  a*  are 
the  values  of  a and  a*  appropriate  for  fully  turbulent  flows 
[see  Equation  (6 ) ] . 

The  value  of  X can  be  fixed  by  demanding  that  the  Saffman- 
model  neutral  stability  Reynolds  number,  Re^^,  for  a Blaslus 
boundary  layer  be  the  sam.e  as  that  predicted  by  linear 
stability  theory.  Neutral  stability  for  the  turbulence 
model  is  defined  as  the  condition  where  a*|9u/9y|  = 6*w 
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^ J where  u Is  turbulent  dissipation  rate^  given  by  u = pS2.  Then, 

using  the  Blasius  velocity  profile  and  a computed  w-profile 
leads  to 

Re^^  - lOOO/X*  (9) 

with  the  least  stable  point  being  at  approj.imately  y/6  = 0.30. 
Hence,  using  the  accepted  linear  stability  value  of  Re^^  of 
9*10®,  X * 0.105.  A similar  argument  for  selecting  the  value 
of  Rq  remains  to  be  found.  Numerical  experimentation  indi- 
cates that  R^  should  be  about  0.10,  but  this  value  should 
o 

be  regarded  as  tentative  until  further  applications  are 
made  to  determine  its  universality. 

Finally,  regarding  the  turbulence  model  and  its  application 
to  transitional  flows,  note  that  algebraic  transition  pre- 
dictions can  be  made  which  apparently  are  accurate  to  within 
a factor  of  two  or  three.  Specifically,  for  many  incompres- 
sible laminar  flows  over  perfectly  smooth  surfaces,  w is 
given  by  0)  « 20v/6y*  so  that  neutral  stability  occurs  when 
y^ I 3u/3y 1 /v=20/(Xa*6/B*)  = 517.  Also,  transition  will  be 
complete  when  production  of  exceeds  w*  dissipation,  i.e., 
when  y^|3u/3y|/v  = 20/(Xa*)  = 722.  Therefore,  an  approximate 
transition  criterion  predicted  by  the  model  for  incompressible 
flows  is 

317  < yM3u/3yl/v  < 722  (10) 

0 

Equation  (10)  resembles  the  Van  Drlest-Blumer  formulation. 


t 

When  dealing  with  incompressible  flows  in  this  study,  the 
turbulent  dissipation  rate,  w,  is  used  in  place  of  the  sped 
flc  turbulent  dissipation  rate,-]. 
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2.3  PREDICTED  TRANSITION  MECHANISM 


Before  discussing  results  obtained  in  the  present  study  for 
constant-pressure  boundary-layer  flows  (Section  3)  and 
blunt-body  flows  (Section  it  is  interesting  to  pause 
and  discuss  the  model-predicted  transition  mechanism.  The 
most-thoroughly  studied  flow  with  regard  to  transition  is 
probably  Incompressible  FPBL  flow;  therefore,  discussion 
below  concentrates  upon  the  Blasius  boundary  layer. 

As  noted  in  Subsection  2.1,  starting  from  the  plate  leading 
edge  with  e = 0 throughout  the  boundary  layer  and  maintain- 
ing e « e^  and  w = at  the  boundary-layer  edge  (y  = 6), 
a small  amount  of  turbulent  energy  is  entrained  beginning 
at  the  plate  leading  edge.  This  turbulent  energy  then 
spreads  through  the  boundary  layer  by  the  action  of  molecu- 
lar diffusion.  At  this  point,  the  turbulent  energy  in- 
creases monatonlcally  from  zero  at  the  plate  surface  to  its 
freestream  value  at  the  boundary-layer  edge. 

No  turbulent-energy  amplification  occurs  for  a significant 
distance  downstream  of  the  plate  leading  edge.  Turbulent- 
energy  entrainment  continues  and.  Just  upstream  of  the 
neutral  stability  point,  a spike,  or  local  maximum.,  deve- 
lops in  the  e-profile  for  y/6  between  0.^5  and  0.60,  de- 
pending upon  the  value  of  This  prediction  is  consis- 

tent with  experimental  m.easurer.ents  which  indicates  that 
disturbances  are  first  amplified  in  a Blasius  boundary 
layer  at  a point  located  about  60^  of  the  way  through  the 
boundary  layer,  ho  effect  on  skin  friction,  shape  factor, 
etc,  is  observed  at  this  point.  The  spike  diffuses  to’.vard 
the  wall  with  the  magnitude  of  the  spike  rem.aining  approx- 
imately constant;  then,  beyond  the  neutral-stability  point, 
the  spike  is  gradually  amplified.  Finally,  at  a nlate- 
length  Reynolds  number  Re^.  , an  abrupt  increase  in  e is 

V 
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observed  throughout  the  boundary  layer,  indicating  the  on- 
set of  transition.  The  spike  is  now  located  at  y/6  » 0.2, 
near  the  critical  layer.  This  prediction  is  consistent 
with  linear  stability  analysis  which  Indicates  the  lease- 
stable  point  in  a Blaslus  boundary  layer  occurs  in  the 
critical  layer.  Turbulent  energy  is  amplified  to  a value 
typical  of  fully  turbulent  flow  and  levels  off  when  Re 

A 

is  two  or  three  times  Rexf  width  of  the  transitional 

region  is  realistic  as  Re^j  is  observed  experimentally  to 
Increase  typically  by  a similar  amount  through  transition. 

Hence,  many  qualitative  features  of  FPBL  transition  are 
adequately  represented  in  the  model-predicted  transition. 
As  will  be  shown  in  the  following  sections,  the  model 
also  yields  accurate  quantitative  predictions  for  proper- 
ties of  engineering  interest  (e.g.,  transition  Reynolds 
number)  in  both  FPBL  and  blunt-body  flows. 
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3.0  CONSTANT-PRESSURE  BOUNDARY-LAYER  FLOWS 


This  study's  primary  objective  has  been  to  develop  an 
accurate  transition-prediction  method  for  blunt-body 
(e.g.,  RV  nosetlps)  flows.  However,  a wealth  of  constant- 
pressure  boundary-layer-transltlon  data  Is  available.  In- 
cluding many  effects  pertinent  to  transition  on  reentry 
vehicles  (RV's)  under  flight  conditions.  Hence,  studying 
both  constant-pressure  boundary-layer  and  blunt-body  flows 
permits  a larger  data  base  for  testing  model  accuracy. 
Additionally,  studying  more  than  one  class  of  flows  may 
permit  development  of  a universally  applicable  transition- 
prediction  method.  Calculations  have  thus  been  performed 
for  PPBL  flow,  channel  flow  and  pipe  flow  in  addition  to 
the  more-pertinent  blunt-body  flows.  Sensitivity  of 
model-predicted  transition  for  PPBL  flow  to  freestream 
turbulence  (Subsection  3.1) » to  surface  roughness  (Sub- 
section 3.2),  and  to  suction  (Subsection  3.3)  Is  deter- 
mined. Channel  and  pipe  flow  are  analyzed  In  Subsection 
3.11. 


3.1  FREESTREAM  TURBULENCE 

With  o and  a*  given  by  Equations  (7  and  8),  effects  of 
freestream  turbulence  on  incompressible  PPBL  transition 
have  been  analyzed.  Table  1 summarizes  the  computations; 
Figure  3 compares  model-predicted  transition  Reynolds 
number,  R^xt*  with  experimental  data^*^^”^^.  Transition 
Reynolds  number  is  defined  in  the  computations  as  the 
point  where  is  first  observed  to  deviate  from  the  lam- 
inar value  by  more  than  0.5!?;  turbulence  intensity,  T*, 

Is  defined  as 


Figure  3*  Transition  Reynolds  nur.ber  for  Inconpresslble 
FPBL  flow  as  a function  of  freestream  turbu- 
lence Intensity,  t' 
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Table  1.  Preestream  turbulence  effects  on  transition 
Reynolds  number  for  Incompressible  PPBL  flow 


As  shown  in  Figure  3,  model-predicted  transition  occurs  at 
much  higher  values  of  Re^  than  were  obtained  when  a*  was 
assumed  Independent  of  Re,j,. 

An  Important  feature  of  the  model  is  that  it  characterizes 
freestream  turbulence  (as  would  be  expected  on  physical 
grounds)  with  two  quantities,  viz,  intensity  and  scale.  In 
addition  to  freestream  turbulent  energy,  the  turbulent  dis- 
sipation rate  w,  must  be  specified;  this  boundary  condition 
is  equivalent  to  specification  of  a turbulence  scale.  The 
effect  of  the  freestream  value  of  w on  predicted  transition 
is  also  shown  in  Figure  3.  The  effect  of  is  most  pro- 
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nounced  for  hlgh-lntenslty  turbulence,  with  Increasing  u 

tending  to  delay  transition.  Excellent  agreement  with  all 

data  shown  Is  obtained  when  vw^/U  ■ .033. 

e e 

3.2  SURFACE  ROUGHNESS 

Surface  roughness  has  been  represented  with  the  turbulence 
17 

model  ' for  fully  developed  turbulent  flows  by  using  the 
following  dissipation  rate  boundary  condition: 

aty«0  (12) 

where  u_  Is  friction  velocity  and  S is  a function  of  surface 

IQ 

roughness.  Saffman  and  Wilcox  have  correlated  S with 
roughness  height,  k,  and  arrived  at  the  relation 

ku^/v  » 50  S“^/^  (13) 

Combining  Equations  (12  and  13)  there  follows 

“ * ^ at  y ■ 0 (1^1) 

00 

Equation  (14)  has  been  developed  with  the  hypothesis  that 
the  roughness  element  Is  very  small  compared  to  a typical 
boundary  layer  dimension  such  as  momentum  thickness,  9.  This 
hypothesis  is  valid  for  roughness  heights  of  practical  inter- 
est (l.e.,  roughness-height  Reynolds  number,  up  to  1000) 

when  a boundary  layer  Is  turbulent.  However,  for  laminar 
boundary  layers,  even  relatively  small  roughness  heights  will 
be  comparable  to  the  boundary-layer  thickness.  For  example, 
when  Rej^  > 150,  a roughness  element  is  more  than  IOTj  of  the 
boundary-layer  thickness  when  Re^  = 10®,  a clear  violation 
of  the  original  postulates  made  in  arriving  at  Equation  (14). 
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A more  general  boundary  condition  has  been  devised  in  the 
present  study,  namely, 

w ■ N(k/0)  ^ (15) 

49 


where 


N(k/e)  1 as  k/0  -*•  0 (l6) 

To  help  deduce  the  functional  dependence  of  N upon  k/0, 
a simple  correlation  between  Rej^  and  the  value  of  u at 
a solid  boundary  has  been  Inferred  from  experimental  data 

7 

of  Pelndt  for  Incompressible  rough-wall  PPBL  flow.  The 
correlation  was  developed  by  numerical  experimentation 
which  showed  that  the  Pelndt  data  are  duplicated  by  model- 
equation  predictions  using  the  generalized  boundary  con- 
dition (15)  with  the  quantity  N given  by 

N * [120/Rej^]*  (17) 

Plgure  ^ shows  the  model  predicted  effect  of  surface  rough- 
ness on  PPBL  transition  when  Equation  (17)  Is  used.  Cal- 
culations were  performed  with  T'  » 1%  and  (vw/U)  ■ .033; 
results  of  the  calculations  are  listed  in  Table  2. 

While  excellent  agreement  with  experimental  data  has  been 
obtained  using  Equation  (17),  the  formulation  is  not  quite 
complete;  that  is,  the  form  of  the  experimental  data  makes 
It  convenient  to  infer  N directly  as  a function  of  Rej^. 
While  Rej^  is  certainly  of  importance  in  rough-wall  tran- 
sition, k/0  Is  a more  fundamental  parame'-er  since  it 
provides  a criterion  for  specifying  when  a roughness 
element  is  large  compared  to  toundary-layer  thickness  while 
the  roughness-height  Reynolds  number  Rej^  does  not.  How- 
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C) 

( R®xt  ^ rough'^  ( R®  xt  ^ smooth 


Figure  4.  Transition  Reynolds  number  for  flat-plate- 

boundary-layer  flow  as  a function  of  Reynolds 
number  based  on  roughness  height,  Re.  . 


ever,  a dependence  of  N upon  k/6  Is  implicitly  contained  In 
Equation  (17)  and  can  be  made  explicit  by  noting  that  for 
laminar  FPBL  flow 


c^Re  » (.66^0* 


wherefore  N can  be  rewritten  as 

I 


N(k/0) 


9 


I < 300c^ 
I < 300Cf 


(18) 


(19) 


Equations  (15  and  19)  are  regarded  as  the  revised  rough- 
wall  boundary  condition  valid  for  both  high-  and  low-Rey- 
nolds-number  boundary  layers.  For  the  compressible  flow 
applications  of  Section  4,  c^  id  defined  In  terms  of  surface 
density,  p , so  that 

W 

Cf  • ZVOw^e  - 


Where  Is  surface  shear  stress. 


Table  2.  Surface  roughness  effects  on  transition  Rey- 
nolds number  for  Incompressible  FPBL  flow 


(Re*t>r 

0 

6.70«10® 

1.00 

108 

6.55*10® 

0.98 

162 

5.95*10® 

0.89 

193 

5.19*10® 

0.77 

219 

4.13«10® 

0.6a 

229 

3.61-10® 

0.54 

261 

1.98*10® 

0.30 

300 

8.25*10- 

0.12 

3.3  SUCTION 


Applications  in  the  preceding  two  subsections  have  served 
mainly  to  fix  the  value  of  the  empirical  constant 
[Equation  (7)]  and  to  determine  the  analytical  dependence 
of  the  function  N upon  k/6  [Equations  (15,  19,  and  20]. 

To  test  the  modified  turbulence  model  with  no  further 
parameter  adjustment,  the  effects  of  suction  on  FPBL 
transition  have  been  analyzed.  Computations  have  been 
performed  for  several  suction  rates  (Table  3)  to  deter- 
mine the  minimum  amount  of  suction  required  to  prevent 
transition;  for  all  calculations  T'  = 2.6^  and  (ufi/U)e 
« .033. 


Table  3.  Effects  of  suction  on  transition  Reynolds 
number  for  Incompressible  FPBL  flow 


Re^ 

^t 

Rev 

^t 

0 

1.0‘10"* 

I.O'IO'^ 

l.i^•10*^ 

9.25‘10'‘ 
1. 05*10® 
1.86-10® 
3.li<-10® 

1.5*10'® 

1.6*10'® 

1.7*10'® 

1.8-10'® 

3.92-10® 

6.03*10® 

> 6.59*10® 

> 6.59*10® 

Results  of  the  calculations  are  shown  in  Figure  5;  Cq 
is  volume  coefficient  defined  as 


Cq  " -V“e  (21) 

where  v^  is  the  suction  velocity.  The  Indicated  minimum 
volume  coefficient  required  to  prevent  transition  is 

<=Qmln  = (22) 

No  transition  occured  in  calculations  having 
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1*10* 


3*10* 


1«10* 


3'10‘* 


.0001 


.001 

(.0017-Cp) 


Figure  5.  Transition  Reynolds  number  for  Incompressible 
flat-plate-boundary-layer  flow  with  uniform 
suction  as  a function  of  volume  coefficient,  Cq. 
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Experimental  data  are  sparse  and  inconclusive  for  transi- 
tional FPBL  flow  with  suction.  Experiments  by  Simpson,  Kays 
and  Moffat^^  for  uniform  suction  indicate  that  lies 

between  .002^1  and  .00^16,  while  Pfenninger *s^®  experiments 
suggest  that  is  within  the  range  .0010  to  .0020.  The 

computed  is  hence  within  experimental  data  scatter.  As 

a final  note,  the  calculated  results  are  far  more  accurate 
than  the  linear  stability  prediction^^  of  = I.IS'KT". 

3.H  CHANNEL  AND  PIPE  FLOW 

Fu?ly  developed  channel  and  pipe  flow  are  especially  simple 
to  analyze  using  the  approximate  transition  criterion  given 
in  Equation  (10).  Table  summarizes  predictions  for  tran- 
sition Reynolds  number,  R,  based  on  average  velocity  and 
channel  height/pipe  diameter;  experimentally  measured  R 
and  the  value  of  R predicted  by  linear  stability  analysis 
are  also  included  in  the  table.  Computed  R agrees  closely 
with  measured  R for  both  flows.  As  with  the  analysis  of 
suction,  model  predictions  are  much  closer  to  corresponding 
measurements  than  are  linear  stability  predictions,  particu- 
larly for  pipe  flow. 


Table  Transition  Reynolds  number  predictions  for 
fully  developed  channel  and  pipe  flow. 


Flow 

R,  Present 

R,  Measured 

R,  Linear 

Analysis 

Stability  Analysis 

Channel 

1427-32^^9 

1400 

7085 

Pipe 

1070-2437 

2300 

00 
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ll.o  BLUNT-BODY  FLOWS 

Results  presented  in  the  preceding  section  for  Incompres- 
sible constant-pressure  boundary-layer  flows  serve  mainly 
to  Improve  the  turbulence  model's  transition  predictive 
accuracy.  In  this  section,  the  model  Is  applied,  with  no 
further  modification,  to  blunt-body  flows.  In  Subsections 
^.1  throgh  ^.3,  effects  of  surface  roughness,  surface 
cooling,  and  body  geometry  are  analyzed  for  sphere-cone 
bodies  Immersed  in  a supersonic  stream.  Subsection  4.^ 
presents  a transition-prediction  correlation  based  on  the 
numerical  results. 

4 . 1 SURFACE  ROUGHNESS 

Two  rounds  of  surface-roughness  calculations  were  performed. 

First,  as  a preliminary  test  of  the  model's  accuracy  for 

blunt-body  flows,  four  cases  from  the  PANT  Series  A V/ind 
2 

Tunnel  Tests  were  simulated.  The  body  considered  in  the 
PANT  Series  A tests  for  the  cases  selected  was  a sphere- 
cone  configuration  with  a nose  radius,  rj^,  of  2.5  in.  Then, 
for  the  same  body  at  flov;  conditions  close  to  those  In  the 
PANT  Series  A experiments,  conditions  for  Incipient  transi- 
tion were  determined. 

4.1.1  PANT  Series  A Computations 

The  cases  considered.  Including  flow  conditions,  are  sum- 
marized in  Table  5.  As  indicated  in  the  table,  roughness 
height,  k,  varies  from  1.5  mils  to  10  mils.  For  all  cases 
the  surface  temperature,  T^,  is  assumed  to  be 

T„/T.  *0.4  (23) 

00 

where  T,.  is  the  freestream  total  temperature.  The  modi- 
^00 

fled  Nev/tonian  pressure  distribution  is  used  to  define 
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boundary- layer-edge  pressure,  l.e.. 


- 1)  + cos*  ♦ (2^») 

where  is  total  pressure  behind  the  shock,  (>  Is  angle 
from  the  centerline,  Is  freestream  density,  and  is 
freestream  velocity.  In  all  calculations,  freestream  tur- 
bulence Intensity  is  T'  ■ 0.33?  while  the  freestream  spe- 
cific dissipation  rate,  0 , is  given  by 

fle  - 0.0033  U*/Mg  (25) 

This  value  for  *2^  is  somewhat  lower  than  the  values  used 

e 

in  the  PPBL  calculations  of  Section  3.  This  lower  value 
is  used  because  larger  values  of  cause,  for  reasons 
unknown,  convergence  problems  near  the  stagnation  point. 


Table  5.  Flow  conditions  for  blunt-body 
roughness  calculations. 


PANT 
Run  No. 

k(mils ) 

00 

Re„(ft"^) 

p^^(psia) 

0 

8 

u 

115 

1.5 

^1.98 

4.10.10^ 

240 

750 

129 

3.0 

**.97 

2.05*10^ 

123 

770 

16H 

10.0 

4.98 

0.55-10^ 

31 

735 

165 

10.0 

4.96 

x> 
0 
1 — ( 

• 

CO 

• 

0 

49 

742 

Transition  location  Is  identified  by  inspecting  computed 
surface  heat  transfer  distribution,  q.  For  example,  Figure  6 
shows  the  computed  variation  of  q for  PANT  Series  A Run  165 
as  a function  of  arc  length, s.  The  ratio  of  q to  its  value  at 
the  stagnation  point,  falls  off  gradually  as  s/r,^ 

increases.  An  abrupt  increase  in  heat  transfer  occurs  near 
s/rjj  = 0.38;  Q/Qgtag  achieves  a maximum  value  of  2.36  at 


s/Pjj  ■ 0.5^»  and  then  gradually  decreases.  For  engineer- 
ing purposes,  the  point  at  which  q first  increases  marks  the 
beginning  of  transition;  similarly,  the  end  of  transition 
could  be  defined  as  the  point  where  q passes  its  maximum 
value. 

Table  6 summarizes  results  of  the  calculations.  Transition 

is  predicted  in  three  of  the  four  cases.  For  the  three 

transitional  cases.  Figure  7 shows  computed  transition 

Reynolds  number  based  on  momentum  thickness,  ^©0^*  a 

function  of  (k/6^)/(T  /T  ) . The  figure  shows  clearly  that 

computed  values  of  Re0  for  these  three  cases  closely  ap- 

t -I  2 

proximate  the  PANT  correlation  * , namely, 

r T *1 

- 215  I (26) 

For  the  one  laminar  run,  PANT  Run  16^,  no  transition  was 
predicted  although  computation  was  carried  out  to  (J>  * ^5^, 
a point  well  beyond  the  sonic  point.  Again,  the  results 
agree  with  measurements,  as  PANT  Run  164  remains  laminar. 

Table  6.  Computed  transition  Reynolds  numbers  for 
blunt-body  roughness  calculations 


PANT 
Run  No. 

k/0^ 

Re« 

®t 

115 

l.i*8 

3.51 

101.7 

129 

2.07 

H.BH 

80.5 

l6i| 

- 

- 

00 

165 

il.50 

10.7^ 

i<4.0 
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4.1.2  Incipient  Transition  Calculations 

Knowledge  of  incipient  transition  conditions  is  important 
in  reentry  vehicle  nosetlp  design.  Hence,  to  provide  a 
further  test  of  the  model’s  ability  to  predict  salient 
features  of  blunt-body  transition,  a series  of  calculations 
was  performed  to  determine,  as  a function  of  roughness 
height,  k,  the  minimum  freestream  unit  Reynolds  number, 

Re^,  at  which  transition  occurs,  l.e.,  the  incipient- 
transition  Re^.  Again,  computations  were  for  flow  near 
the  stagnation  point  of  a sphere-cone  body  having  a nose 
radius,  r^j,  of  2.5  inches.  For  simplicity,  in  all  cal- 
culations freestream  Mach  number,  M^,  and  freestream  total 
temperature,  T.  , were  5 and  750®F.,  respectively;  surface 
temperature  was  0.4t.  . The  values  of  T’  and  (yfi/U*) 

6 

were  0.3358  and  0.0033. 

Pour  roughness  heights  were  considered,  namely,  0.5  mil, 

1.5  mils,  3 mils  and  10  mils.  Table  7 summarizes  results, 
including  momentum  thickness  at  the  stagnation  point,  Qg^-^gi 
transition  point  values  of  angle,  (J)^,  boundary-layer-edge 
Mach  number,  M , displacement  thickness,  6f,  and  momentum 
thickness,  6.  ; and  the  PANT  coordinates,  (k/9. )/(T  ,/T^) 

Xt  u W 6 

and  Re„  . 

®t 

Figure  8 illustrates  calculated  behavior  for  k » 3 nils. 

Each  symbol  denotes  transition-point  location.  As  shown, 
four  of  the  points  lie  along  the  numerical  transition  cor- 
relation (NTC)  line  defined  by 


which  is  very  close  to  the  PANT  correlation  [Equation  (26)]. 
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■"i  Table  7»  Summary  of  Incipient-transition  surface-roughness  calculations. 

1 } 
s y 


— 

6*(mils) 

0^(mils) 

k / 

»t  Te 

Re  a 
®t 

0 

.5 

8.il2*10^ 

n 

21.5 

0.48 

M 

0.70 

1.71 

134 

6.00*10° 

27.2 

0.61 

0.86 

1.35 

142 

4.00*10° 

0.96 

40.2 

0.96 

0.84 

1.20 

0.88 

164 

3.70*10° 

1.00 

46.2 

1.15 

1.05 

1.36 

0.73 

177 

3.60*10° 

1.01 

49.8 

1.27 

1.20 

1.46 

0.65 

186 

3.55-10° 

1.02 

52.8 

1.39 

1.35 

1.57 

0.57 

193 

3.50-10° 

1.02 

- 

- 

- 

- 

- 

00 

1. 

5 

4.00*10^ 

0.96 

23.9 

0.53 

0.60 

1.03 

3.45 

103 

3.00*10° 

1.11 

36.3 

0.85 

0.88 

1.32 

2.48 

130 

2.70*10° 

1.17 

51.7 

1.34 

1.49 

1.76 

1.56 

166 

2.65-10° 

1.18 

- 

- 

- 

- 

- 

00 

3. 

0 

2.00*10^ 

1.36 

27.8 

0.63 

0.91 

1.50 

4.64 

83 

1.75-10° 

1.45 

35.2 

0.82 

1.13 

1.70 

3.88 

97 

1.65-10° 

1.49 

42.6 

1.03 

1.42 

1.92 

3.22 

111 

1.63-10° 

1.50 

47.3 

1.18 

1.64 

2.09 

2.81 

120 

1.62-10° 

1.50 

56.4 

1.54 

2.34 

2.56 

1.99 

137 

1.60-10° 

1.51 

- 

- 

- 

- 

- 

CO 

10. 

0 

0.80-10^ 

2.15 

23.4 

0.52 

1.34 

2.29 

10.34 

45 

0.70-10° 

2.28 

27.8 

0.63 

1.55 

2.53 

9.15 

49 

0.60-10° 

2.46 

36.0 

0.84 

1.98 

2.94 

7.45 

58 

0.58*10° 

2.52 

39.6 

0.94 

2.20 

3.12 

6.49 

62 

0.55*10° 

2.59 

- 

- 

- 

- 

- 

00 
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Illustration  of  the  method  for  locating  incipient 
transition  freestream  unit  Reynolds  number; 

;i  = 3 mils;  = O.li. 


The  point  for  Re^  » 1.62»10^ft”^  Is  the  incipient  case 
and  lies  slightly  below  the  NTC  line.  The  three  curves 
(referred  to  as  Re^  trajectories)  are  based  on  computed 
values  of  Re-  and  (k/0)/(T  /T  ) up  to  transition.  On  the 
one  hand,  for  Re^  « 2*10°ft“  the  Re^  trajectory  inter- 
sects the  NTC  line,  and  transition  occurs  at  the  inter- 
section; on  the  other  hand,  when  Re^  = 1.60*10^ft“^,  the 
Re^  trajectory  never  Intersects  the  line  and  no  transi- 
tion occurs.  The  incipient  case  lies  between  these  two 
extremes.  As  shown  in  Figure  8,  the  Re^  trajectory  is 
tangent  to  the  NTC  line.  However,  transition  does  not 
occur  at  the  point  which  would  be  about  midway  between 
the  NTC  line  and  the  PANT  correlation  curve.  This  be- 
havior is  consistent  with  claims  that  transition  loca- 
tion and  transition  onset  (i.e.,  incipient  transition 
location)  are  not  coincident. 

Figure  9 shows  that,  .•;ith  the  exception  of  the  0.5-mil 
case,  all  the  computed  transition-point  locations  listed 
in  Table  7 lie  along  the  NTC  line.  Since  corresponding 
experimental  data  also  defy  correlation,  the  fact  that 
0.5-mil  numerical  data  fail  to  correlate  lends  further 
confidence  to  the  predictions.  As  with  the  3-mil  case, 
the  incipient  transition  point  lies  below  the  NTC  line 
when  k = 1.5  mils  and  k = 10  mils.  Figure  10  presents 

computed  incipient-transition  Re_  as  a function  of  k; 

2 

computed  values  are  very  close  to  measured  values  . 
l\.2  SURFACE  COOLING 

As  with  surface  roughness,  effects  of  surface  cooling  on 
transition  for  the  r,^  = 2.5-in.,  sphere-cone  body  have 
been  analyzed  in  two  steps.  In  the  first  step,  k and  Re^ 
remain  constant,  and  is  varied.  In  the  second  step, 
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.3  1 3 10 

(k/e^^)/(VTe) 


Figure  9.  Summary  of  surface  roughness  and  surface  temperature 
calculations  for  sphere-cone  stagnation-region  flow. 
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Figure  10. 


Freestream  unit  Reynolds  number  for  incipient 
transition  as  a function  of  roughness  height; 


T„/T, 


= OA 


with  Re_  held  constant,  k is  varied  for  T„/T.  « 0.2 

• W too 

and  0.8  until  Incipient  transition  conditions  are  deter- 
mined. 


For  the  first  set  of  calculations,  the  roughness  height 
Is  3 mils  and  the  corresponding  incipient-transition  Re^ 
for  T /T^.  = 0.^1  of  1.62«10®  Is  used.  As  summarized  in 
Table  8,  surface  temperature  Is  varied  from  0.1  to  0.8; 
computed  Re 9^  is  shown  graphically  in  Figure  11  as  a func- 
tion of  T /T^  . As  would  be  consistent  with  experimental 
measurements  for  sphere-cones,  surface  cooling  is  predict- 
ed to  have  a destabilizing  effect  on  transition.  Figure  9 
shows  that  for  T, ,/T^  > 0.2,  predicted  transition  location 

Is  close  to  the  PANT  correlation  line.  However,  for 
T,,/T4.  = 0.1,  ROfl^  is  about  ^10^  lower  than  the  PAHT-cor- 

relation  value. 


In  the  second  set  of  surface-cooling  calculations,  Re^ 
is  1.62*10*;  k is  then  varied  for  T /T^  = 0.2  and  0.8 

until  incipient-transition  conditions  have  been  determined. 
Results  of  the  calculations  are  summarized  in  Table  8 
and  Figure  12.  An  approximate  fit  to  the  numerical  data 
for  the  Incipient-transition  roughness  height,  is 


‘‘Incip 


10  (VT 


(28) 


with  given  in  mils 


. 3 BODY  GEOMETRY 


Effects  of  body  geometry  were  analyzed  by  computing  in- 
cipient-transition conditions  for  sphere-cones  of  varying 
nose  radius.  For  all  of  the  calculations,  the  following 
conditions  have  been  used; 

= 5 '^t  ” I 

T^  = 500°  R.  k ” = 3.5  mils  | (29) 

T’  = 0.33s  = 0.0033  Ug/Ug  I 
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Table  8,  Summary  of  surface-cooling  calculations. 


k(mlls ) 
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6*(mils) 

6^(mils) 
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0.10 

3. 

00 

1.67 

6.6 

0.14 

0.42 

1.68 
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21 

0.20 

1.62 

13.8 

0.30 

0.03 

1.67 

8.90 

41 

0.30 

1.56 

23.4 

0.52 

0.51 

1.66 

5.69 

66 

0.38 

1.51 

36.6 

0.86 

1.11 

1.81 

3.80 

97 

0.45 

1.48 

- 

- 

- 

- 

- 

00 

0.80 

1.30 

- 

- 

- 

- 

- 

00 

0.20 

2. 

00 

1.62 

20.4 

0.45 

0.06 

1.70 

5.65 

60 

1. 

50 

28.1 

0.63 

0.10 

1.80 

3.86 

82 

1. 

00 

50.3 

1.29 

0.45 

2.44 

1.54 

138 

0. 

98 

53.6 

1.42 

0.55 

2.63 

1.33 

146 

0. 

97 

56.1 

1.53 

0.65 

2.81 

1.18 

152 

0. 

95 

- 

- 

- 

- 

00 

0.80 

10. 

00 

1.30 

27.0 

0.61 

2.55 

1.40 

8.29 

62 

8. 

00 

36.9 

0.86 

2.99 

1.51 

5.76 

82 

7. 

50 

52.3 

1.37 

4.49 

1.86 

3.67 

105 

7. 

49 

- 

- 

- 

- 

- 

00 

Boundary-layer-edge  pressure  is  again  given  by  the  modified 

Newtonian  distribution  [Equation  (2^)].  Three  nose  radii  are 

considered,  namely,  0.75  in.,  2.5  in.  and  3-5  in.;  Table  9 

summarises  the  calculations.  As  shown  in  Figure  13,  rj^  has 

only  a slight  effect  on  the  incipient-transition  with 

Increasing  nose  radius  yielding  lower  values  for  Re„.  This 

31 

trend  is  consistent  with  recent  measurements  , thus  providing 
another  test  of  the  theory.  The  numerical  data  are  closely 
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k(mils) 

Figure  12.  Minimum  roughness  height  for  incipient  transition 

as  a function  of  surface  temperature;  Re^  = 1.62*10 
per  foot. 
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approximated  by  the  following  formula: 

Re„  « 1.48*10^r  ft"^  (30) 

00  ig 

with  r^^  given  In  Inches. 

il.il  NUMERICAL  NOSETIP-TRANSITION  CORRELATION 


All  computed  Incipient  transition  points  are  shown  In 
Figure  li|.  Six  of  the  nine  points  correlate  with  the 
PANT  curve  [Equation  (26)].  As  In  the  PANT  experiments, 
the  0.5-mll  data  point  Is  one  of  the  points  which  Is  not 
represented  adequately  by  the  curve.  The  other  two  points 
fall  to  correlate  because  the  predicted  effect  of  surface 
temperature  Is  stronger  than  that  Indicated  by  the  PANT 
curve.  Figure  15  presents  a new  correlation  which  works 
well  for  all  but  the  0.5-mll  numerical  data  point.  Spe- 
cifically, the  new  correlation  Is 


(31) 


-m- 


Figure  1^1.  Comparison  of  computed  Incipient  transition 
points  with  the  PANT  correlation. 
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Figure  15.  Incipient  transition  correlation  inferred  fror. 

numerical  results  for  sphere-cone  stagnation- 
region  flow. 
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5.0  DISCUSSION 


Results  presented  in  Sections  3 and  ^ show  that  accurate 
transition  predictions  can  be  made  with  the  Saffman- 
turbulence-model  equations,  provided  low-Reynolds-number 
effects  are  taken  into  account  (Section  2).  V/ith  only 
one  adjustable  parameter,  namely  R^,  the  model  yields  ac- 
curate transition  predictions  for  incompressible  boundary 
layers  and  for  stagnation-region  flows.  The  effects  of 
freestream  turbulence,  surface  roughness,  and  suction  on 
Incompressible  FPBL  transition  are  well  represented  by 
the  model;  surface-roughness,  surface-cooling,  and  body- 
geometry  effects  on  blunt-body  transition  also  have  been 
predicted  accurately. 

The  fact  that  the  model  works  equally  well  for  constant- 
pressure  boundary  layers  and  for  stagnation-region  flows 
is  a key  feature  of  the  technique.  That  is,  although 
further  study  will  be  needed  to  establish  the  model’s 
range  of  applicability,  the  fact  that  the  initial  data 
base  includes  these  two  classes  of  flows  provides  hope 
that  a universally  applicable  transition-prediction 
method  has  been  developed. 

Specific  areas  needing  further  model  development  and 
validation  include  effects  on  transition  of  freestream 
turbulence,  compressibility  and  pressure  gradient.  As 
shown  in  Section  2,  freestream-turbulence  scale  has  a 
large  effect  on  FPBL  transition.  V/hile  some  effect  is 
expected  on  physical  grounds,  no  attempt  has  been  made 
to  determine  whether  or  not  the  predicted  effect  is  even 
qualitatively  correct.  Sensitivity  of  blunt-body  tran- 
sition to  both  freestream-turbulence  intensity  and  scale 
should  also  be  determined.  The  predicted  role  of  free- 
stream turbulence  in  rough-wall  transition  requires 
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further  elucidation.  Compressibility  effects  should  be 
analyzed  In  detail  for  FPBL  flows  for  which  extensive 
experimental  data  are  available.  Surface  cooling  and 
freestream  Mach  number  are  easily  represented  with  the 
model  and  should  be  analyzed  to  help  test  validity  of 
the  predicted  stronger-than-measured  effect  of  surface 
cooling  on  blunt-body  transition.  While  pressure-gradient 
effects  have  been  addressed  Indirectly  by  analyzing  FPBL 
flow  and  the  blunt-body  geometric  configurations,  more- 
direct  tests  of  the  model  are  needed.  Again,  extensive 
experimental  data  exist  for  the  effect  of  pressure  gra- 
dient on  boundary-layer  transition,  so  that  definitive 
tests  of  the  model  can  be  easily  made. 

In  conclusion,  the  computations  In  Section  K exemplify 
the  turbulence-model-transltlon-predlctlon  method's  poten- 
tial value  to  the  reentry-vehlcle-nosetlp-deslgn  engineer. 
The  EDDYBL  computer  code  has  been  used  in  Section  as  a 
"numerical  wind  tunnel"  In  that  PANT  Series  A experiments 
have  been  simulated  Including  a wider  range  of  surface 
temperatures  than  could  be  covered  in  the  physical  tunnel. 
Many  of  the  predicted  results  agree  with  measured  flow 
properties.  In  particular,  measured  and  computed  effects 
of  surface  temperature  on  sphere-cone  transition  are  In 
close  agreement  for  the  range  of  temperatures  in  the  PANT 
experiments;  however,  computations  at  lower  temperatures 
than  those  which  could  be  considered  in  the  PANT  experi- 
ments indicate  that  surface  cooling  has  a stronger  effect 
than  that  Inferred  from  the  experimental  data.  This  in- 
dicates an  area  in  which  further  measurements  are  needed. 
Additionally,  this  exemplifies  the  fact  that  arbitrary 
flow  conditions  can  be  simulated  In  the  numerical  wind 
tunnel  while  the  physical  wind  tunnel  often  Is  limited. 
Hence,  the  numerical  wind  tunnel  can  be  used  to  test  and 
verify  the  transition-prediction  method  for  both  ground- 
test  and  flight-test  data  and,  as  the  ultimate  goal,  can 
be  used  to  predict  nosetlp  transition  under  real  flight 
conditions . 
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APPENDIX 

THE  EDDYBL  COMPUTER  CODE 

Computations  in  this  study  were  performed  with  the  EDDYBL 

computer  code.  This  code  is  based  on  a boundary-layer 

21 

program  developed  at  the  NASA  Langley  Research  Center  ; 

the  Saffman  turbulence  equations  were  incorporated  in  the 

30 

code  in  a DCW  Industries  research  project  . The  integra- 
tion method  embodied  in  EDDYBL  is  similar  to  the  implicit 

31 

Plugge-Lotz  and  Blottner  technique  in  which  the  momentum 
and  energy  equations  are  coupled.  In  adding  the  turbulence 
model  to  the  code,  the  two  nonlinear  diffusion  equations 
(^,  5)  are  solved  in  a coupled  manner  analogous  to  that 
used  for  the  momentum  and  energy  equations  in  the  origi- 
nal version  of  the  code.  Since  all  four  equations  are 
not  solved  simultaneously,  an  iterative  procedure  is 
needed  to  obtain  an  accurate  solution.  As  part  of  this 
study,  a more  efficient  iterative  procedure  than  had  been 
used  previously  was  developed. 

The  new  iterative  scheme  (Figure  Al)  takes  advantage  of 
the  special  nature  of  transition-prediction  calculations. 
Specifically,  note  that  for  laminar  flow  the  momentum  and 
energy  equations  are  not  coupled  to  the  turbulence-density 
equations.  Hence,  with  no  sacrifice  of  accuracy,  the 
momentum  and  energy  equations  can  be  solved  non' teratively 
at  each  station  on  laminar  regions.  Since  the  turbulence 
density  equations  do  require  iteration  when  the  flow  is 
laminar,  eliminating  the  momentum  and  energy  equations 
from  the  iterative  loop  significantly  reduces  computing 
time.  In  specifying  when  all  four  equations  must  be 
solved  3 tci'atively , the  criterion  is  that  the  maximum 
value  (with  respect  to  distance  from  the  boundary)  of  the 
kinematic  eddy  viscosity,  e,  exceed  5^^  of  the  kinematic 


Figure  Al. 


Schematic  of  iterative  method  used  in  the  EDDYBL  code 


molecular  viscosity,  v. 


Experience  with  EDDYBL  shows  that  the  turbulent  energy 
converges  more  slowly  than  the  other  field  variables. 

Two  of  the  most  sensitive  turbulent-energy  values  are  the 
peak  turbulent  energy  and  the  value  at  the  mesh  point 
nearest  the  surface.  These  two  quantities  are  therefore 
monitored  to  determine  when  convergence  Is  attained; 
convergence  Is  defined  to  occur  when  these  two  quantities 
change  by  less  than  \%  between  successive  Iterations. 

Two  or  three  Iterations  appear  to  be  suitable  on  both 
laminar  and  turbulent  regions,  while  seven  or  eight 
Iterations  typically  are  required  through  transition. 
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